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The interplay between disorder and superconductivity has intrigued physicists 
for decades. Of particular interest is the influence of disorder on the supercon- 
ducting energy gap A. In the absence of Coulomb interactions between electrons, 
disorder leads to emergent granularity of the local order parameter resulting in 
a pseudogap at temperatures above the critical temperature T c , as well as a fi- 
nite gap A on the insulating side of the disorder-driven superconductor-insulator 
transition (SIT). At the same time, disorder also enhances the Coulomb inter- 
actions, which subsequently may influence A in a manner that is still not fully 
understood. Here we investigate the evolution of the energy gap through the 
SIT by two different experimental methods: tunneling spectroscopy, in which 
a metallic electrode is placed close to the studied sample thus screening the 
Coulomb interactions, and terahertz (THz) spectroscopy, which probes the un- 
screened sample. The comparison between the two methods illustrates the role 
played by electronic interactions in determining the nature of the phases across 
the SIT and sheds light on the mechanisms involved in the destruction of super- 
conductivity. 

Introduction. The effect of disorder on superconductivity, even for s-wave superconduc- 
tors, is an outstanding problem, yet unresolved. It involves the interplay of three major 
phenomena: superconductivity describing the pairing of electrons due to retarded attractive 
interactions, localization of wave functions because of strong disorder, and the evolution of 
Coulomb interactions with increasing disorder. It was argued by Anderson 1[ that supercon- 
ductivity is not affected by disorder on the scale of the superconducting energy gap h/r « A, 
where r _1 is the scattering rate, persisting even in polycrystalline or amorphous materials. 
Experimentally it was found that superconductivity can be destroyed by a sufficiently large 
degree of disorder of the order of the Fermi energy 2ni| h/r ~ €f. Once superconductivity 
is destroyed the sample undergoes a transition to an insulating state across a superconductor 
insulator transition (SIT), a fundamental manifestation of a quantum phase transition at 
T = tuned by a non-thermal parameter such as disorder or a magnetic field. 

In this paper we focus on the effect of disorder on the energy gap A. According to the 
BCS theory the relationship between the critical temperature T c and the energy gap A is 



given by 2A//cbT c = a where the weak coupling BCS result gives a ~ 3.5. The Mcmillan 
formula [8| in the strong electron-phonon coupling limit provides no such explicit formula but 
experimentally a is found to increase with increasing coupling strength reaching a typical 
value of 4.8 for strong coupling superconductors such as Pb 9|. 

Disorder may influence the energy gap in various ways. Finkelstein [l0(] showed that 
in uniform superconductors, disorder induces diffusive transport which in turn intensifies 
Coulomb repulsion between electrons. This suppresses electron pairing leading to a decrease 
of both the critical temperature and the energy gap with increasing disorder. Indeed a thin 
normal metal was shown to increase T c and A of a nearby disordered superconductor by 
screening the coulomb interactions [ill ]. 

Another possibility is that the superconductor, which exhibits a many-body coherent 
state of pairs on the scale of the coherence length £ 3> kp 1 (fop is the Fermi wave vector) 



becomes spatially inhomogeneous 12j with increasing disorder. This gives rise to an emer- 



gent granularity of the pairing am plit ude on the scale of £ even when the random potential 



fluctuates on an atomic scale 



13Ml5J]. It has been further shown that phase fluctuations 



riding on top of the inhomogeneous pairing amplitude suppress the superfluid stiffness D 



at a critical disorder, signaling the superconductor to insulator transition [161 ] . However, the 
energy gap remains robust and finite across the transition and also survives as a pseudogap 
above the critical temperature T c . Indeed, recent tunneling measurements on InO and TiN 
films have confirmed these theoretical predictions and showed that as disorder increases and 
the sample approaches the SIT, a increases monotonically and ultimately even diverges at 



the SIT where T c — >• but the energy gap A is sti 
gap was observed above the critical temperature 



18 



1 finite. Consequently, a finite pseudo- 



disorder in the insulating phase 21] 



201 ] and a hard gap beyond the critical 



An important question in the field remains open: how does the combination of both 
disorder induced "granularity" and strong Coulomb interactions affect the energy gap? Does 
the joint contribution lead to a constant a or a change of a with disorder? We employ 
two different experimental methods to probe disordered superconductors and compare the 
results of the energy gap: The first is tunneling spectroscopy which requires a close metallic 
electrode that naturally induces screening and thus suppresses Coulomb interactions. The 
second is THz spectroscopy which utilizes free space radiation through the sample, with no 
metallic electrodes or contacts necessary thus enabling the investigation of the unscreened 



unperturbed sample. As a cross-check we also perform THz spectroscopy measurements 
in the presence of a metal electrode. Our main observations are: (1) When the Coulomb 
interactions are screened there is a finite gap in both the superconductor and the insulator, 
albeit with reduced coherence peaks in the insulator. (2) When the Coulomb interactions 
are unscreened the gap in the disordered superconductor is substantially suppressed and in 
the insulator no discernible gap is detected. Furthermore, the gap value A is suppressed 
with disorder much faster than T c and consequently a decreases monotonically to as the 
sample is driven through the SIT. 

Paradigms for the SIT. The disorder driven SIT can be conceptualized in terms of three 
paradigms depicted in Fig. [TJ These are based on the behaviour of energy scales obtained 
from the interplay of BCS superconductivity, Anderson localization and Coulomb correla- 
tions. Paradigm (a) shows the behaviour of the energy gap A and the superfluid stiffness 
D s with increasing disorder, based on calculations that include inhomogeneous pairing am- 



plitude and phase fluctuations [13|, [14], [16|] but ignore Coulomb correlations between the 
electrons. D s = (h 2 /m)n s (where m is the Cooper pair mass and n s is the superfluid den- 
sity) has dimensions of energy for two-dimensional superconductors and can be interpreted 
as the energy associated with phase fluctuations. T c is determined by the minimum of the 
two energy scales. The vanishing of D s at a critical disorder signals a quantum phase transi- 
tion from a superconductor to a gapped Bose insulator, where Cooper pairing still survives 



in the insulator in the form of localized Cooper pairs [17| , as also confirmed by tunneling 



experiments 181421]. 



la. 



Panels (b) and (c) illustrate the additional effects of Coulomb interactions and screening, 
or rather the lack of it, with increasing disorder. Depending on the interplay between A and 
D s there are two possibilities that govern the nature of the phases and the transition. In 
paradigm (b), there is a large regime where A in the presence of both disorder and Coulomb 
interactions is finite. Once D s dips below A, the transition temperature T c is determined 
by D s , the minimum of the two scales. This results in a transition from a superconductor 
to a Bose insulator driven by vanishing D s and a further crossover from a Bose-insulator of 
localized Cooper pairs to a Fermi insulator of broken pairs. This paradigm may be relevant 



for the observed giant magnetoresistance peak in several disordered superconductors [22H30]. 
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FIG. 1: Paradigms for the disorder driven superconductor-insulator transition (SIT). 
Schematic description of three paradigms (see text) for the effect of disorder on the characteristic 
energy scales of a superconductor: the energy gap A (red line) and the superfluid stiffness D s 
(dark blue line). The interplay between these energies defines three phases: a superconductor (SC) 
containing free cooper pairs, a Bose insulator (BI) containing localized pairs and a Fermi insulator 
(FI) containing localized quasi-particles. 

In paradigm (c) enhanced Coulomb interactions suppress A to zero and consequently also 
D s which vanishes at the same critical disorder. In this case one can expect a transition 
from a superconductor to a Fermi insulator directly; this indeed appears to be the case in 
the present THz results presented below. 

We next discuss our experimental results and return to a discussion of what they imply 
in the context of the above paradigms. 



Experimental results and discussion. We have employed tunneling spectroscopy and 
THz optical measurements as the prime methods to determine the energy gap of disordered 
superconductors. To that end, the tunneling spectra obtained at certain temperatures are 
fitted by the BCS density of states; the complex THz conductivity is described by the BCS- 
based Mattis-Bardeen theory [3ll. l32| (see methods section). Superconductivity is reflected 
in the characteristic behaviour of the complex conductivity (cri+i^) as a function of photon 
energy hf, where / the frequency and h Planck's constant. For frequencies / = 2A/h a kink 
in <Ji(f) is observed due to the onset of radiation absorption; a 2 diverges as 1/ / indicating 
the presence of a finite superfluid stiffness D s . Hence this method has the additional value 



that it allows the extraction of both A and D s . We have utilized these two techniques to 
determine the behaviour of A through the SIT in two disordered superconducting materials: 
amorphous indium oxide (InO) films and niobium nitride (NbN) films. In both materials the 
chemical composition can be manipulated thus enabling preparation of batches of samples 
with different disorder strengths, characterized by the normal state sheet resistance R n of 
the films spanning the SIT. In both systems a pseudogap has previously been observed in 
the tunneling DOS for T > T c 19|, |20|. In InO a similar gap was also observed for the 
insulating phase of the SIT [21 ]. 
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FIG. 2: THz versus tunneling spectroscopy on insulating and superconducting indium 
oxide films. T c = 3.3K. a, The real part of the conductivity (upper panel) and the imaginary 
part of the conductivity (lower panel) obtained by using two different methods: BCS model (solid 
curves) and a manual calculation based on Fresnel equations (dots with error bars), see methods 
section for more details. Note that while the superconducting film exhibits a Mattis-Bardeen 
behaviour typical for superconductors, the insulating film exhibits constant behaviour typical for 
normal state, b, Tunneling results on similar indium oxide films. Both the insulating and the 
superconducting films exhibit a similar energy gap [2l[ |. the surrounding depleting background 



indicates the degree of localization 



33| 



Figure |2] compares the THz and tunneling results for two InO samples close to the SIT, one 
on the superconducting side and the other on the insulating side. The THz measurements 
reveal a finite gap (2A = 0.62 meV) and finite superfluid stiffness (D s = 7 meV) in the 
superconducting phase but there is no sign of a gap nor for superfluid stiffness down to 
~ 0.05 meV for the insulator. On the other hand, tunneling measurements conducted on 
the same two samples exhibit a finite gap on both sides of the transition. Though neither 
curves fit the BCS tunneling DOS very well, the best fits yield A « 0.7 meV for both the 
superconductor and the insulator (see ref. |2l|). Similar trends were observed in previous 
planar junctions and STM measurements on both InO and NbN films 19L l20j . Evidently, the 
presence of a tunneling electrode not only enhances the gap in the superconducting phase, 
but in addition, the gap remains finite in the insulating phase. 

This trend is further demonstrated in Fig. [3] which depicts 2A/A;bT c obtained by tun- 
neling and by THz spectroscopy in batches of NbN and InO films having different degrees 
of disorder. For relatively clean samples, both methods yield the same A and correspond- 
ingly the same a = 2A/ksT c . As disorder becomes stronger, T c decreases and eventually 
approaches zero, however, tunneling measurements always observe a finite gap, resulting 
in an increase of a with disorder. The results obtained from THz spectroscopy are quite 
the opposite. With increasing disorder, an energy gap can still be extracted from the con- 
ductivity data yielding good fits to BCS theory; however its energy scale is significantly 
suppressed compared to the BCS expectations. In both InO and NbN films a decreases 
rapidly with disorder in stark contrast to the trend observed in tunneling experiments. In 
fact, A becomes unmeasurable in the THz frequency range well before T c reaches 0. 

The fact that the THz and tunneling spectroscopy agree very well for clean films but 
show opposite trends in the strongly disordered limit demonstrates the crucial role played 
by Coulomb interactions. The tunneling electrode influences the measurement by screening 
long range interactions, thus probing a sample in which interactions are fully or partially 
"switched off". In this limit a increases with disorder. The THz experiment measures an 
unperturbed sample with Coulomb interactions present. The "unscreened" Athz is clearly 
smaller than the "screened" A tunne i ing and the difference between the two becomes larger 
with increasing disorder. 

The inset of Fig. [3] depicts the two energies A and D s , extracted from the THz data, as 
a function of disorder represented by R n - Both scales drop sharply to zero with increasing 
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FIG. 3: THz versus tunneling spectroscopy for series of NbN and InO films with 
different disorder, a = 2A/ksT c as a function of sheet resistance for batches of InO films (left) 
and NbN films (right) obtained by planar junction tunneling (red squares), STM tunneling (red 
circles), and THz spectroscopy (blue circles). The dashed line indicates the experimental value of 



lg] and 



2CJ]. Note that the results 



clean samples. Some of the tunneling results were taken from 
from planar junction and STM methods agree with each-other very well. Inset: The energy gap 
and the superfluid stiffness as a function of sheet resistance for the NbN batch. The resemblance 
to paradigm c (Fig. [lb) is evident. 

disorder in a way that clearly resembles Fig. [Tfc. This behavior of the gap and the superfluid 
stiffness support the notion that the samples experience a transition from a superconductor 
directly to a Fermi insulator. 

Next we present results from the cross-check experiment in which the gap is determined 
from the THz experiment on films with and without a nearby metal plate. For this purpose 
we fabricate two InO films which are deposited together under identical conditions, one is a 



"bare" sample on an insulating substrate and the other a "screened" sample separated from 
an Al metal layer by a 5 nm thick insulator. As seen in the inset of Fig. HI the presence 
of the metal plate increases the T c of the film. a-i measured at T c bare < T < j' c screened for 
both samples is shown in Fig. HI While the bare sample shows no change in cr 2 (f) and 
remains constant at these energies as expected for the normal state, the screened sample 
exhibits a divergence of <Ji at low / indicating the presence of a non-zero D s and a finite 
superconducting gap A ~ 0.1 meV. Note that in this setup the direct measurement of <7\ 
is less accurate due to the absorption of the metal film which drastically cuts down the 
transmission; nevertheless, the phase shift and consequently o<i are measured reliably. From 
the above results we make two important conclusions: First, unlike tunneling results, THz 
measurements show no gap for T > T c and secondly, at T = 1.85K the presence of the 
nearby metallic plane induces a gap in the film that is otherwise insulating and gapless, thus 
revealing the role of Coulomb interactions in suppressing A. 

Still, the fact that in the presence of interactions disorder affects A much more than T c is 
puzzling. According to Fig. [T^c) T c should be determined by A, the smaller of the two energy 
scales. But why does a decrease significantly with disorder? The prediction of Finkelstein 
[lOj captures the suppression of the pairing gap due to enhanced Coulomb interactions. 
However, it does not explain why large disorder suppresses the gap much stronger than T c . 
While a detailed theory still needs to be fleshed out one may consider the influence of the 
normal state DOS. Coulomb interactions play a significant role in restructuring the DOS 
of a disordered normal metal close to the metal-insulator transition. For weakly localized 
metallic films the interplay between Coulomb interactions and disorder leads to a depletion 
of the DOS around the Fermi level 33j in what is known as the zero bias anomaly. For 
strongly localized systems, electron interactions open a soft Coulomb gap at the Fermi energy 
as shown by Efros and Shklovskii [34]. Indeed, all tunneling measurements on disordered 
superconductors exhibit a normal-state background showing suppression of the density of 



state around zero energy |20j, |21]. Hence, the superconducting energy gap must open up 
on this background of a strongly depleted density of state. It is reasonable that this may 
have a considerable effect on the superconducting gap energy. The THz experiments imply 
that this is indeed the case and we expect them to invoke theoretical ideas on the interplay 
between disorder and interactions that would shed light on the observed results. 
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FIG. 4: THz spectroscopy of an indium oxide samples with and without a nearby 
metal plane. The imaginary part of the conductivity for the bare sample (black curve) and 
the screened sample (dark yellow curve) obtained at T = 1.85-ftT and normalized to the results at 
T = 2.5K. The rapid increase in the imaginary part, seen only in the screened sample, indicates the 
superconducting DC conductivity by virtue of Kramers-Kronig relation. The bare sample shows 
constant behaviour as expected from the normal state. The difference between the two identical 
InO samples is further illustrated in the inset which shows the resistance versus temperature for 
both samples. The screened sample exhibits a higher critical temperature than the bare sample. 

Conclusions The set of results presented here reveal the important role Coulomb cor- 
relations play on the energy gap of disordered superconductors. In the absence of Coulomb 
interactions, disorder suppresses D s and hence T c much more than the gap. This is in- 
deed what is observed by the tunneling measurements (in which interactions are naturally 
screened due to the presence of a metallic electrode or a tip) where a gap is observed for 
T > T c as well as in the insulator. Thus the tunneling measurements are well described 
by paradigm (a). Unscreened Coulomb interactions produce the opposite trend as revealed 
by THz results. The surprising finding is that a = 2A/KbT c decreases monotonically with 
increasing disorder indicating that A is suppressed by disorder more than T c . This can not 
be accounted for by the existing paradigms for the SIT. The results also show that when 
Coulomb interactions are screened by placing a metal electrode a gap is recovered in the 
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insulator. This is a very significant result suggesting the possibility of tuning a quantum 
phase transition from a gapless Fermi insulator to a Bose insulator by enhanced screening 
of the Coulomb interactions. Further theoretical studies that incorporate the effects of both 
the emergent granularity and Coulomb interactions are essential for explaining our exper- 
imental results and for a better understanding of when the different paradigms discussed 
above are borne out in experiments. 



Methods 

The InO films were deposited on 10x1 0mm 2 of THz-transparent MgO or sapphire sub- 
strates (with various thickness ranging from 0.5 to 1.5 mm) by e-gun evaporation. During 
the deposition process dry oxygen was injected into the chamber; the partial oxygen pressure 
allows us to tune the disorder. The NbN films were grown on similar MgO substrates by 
reactive magnetron sputtering, where the Nb/N ratio in the plasma served as a disorder 
tuning parameter. In both cases the deposited films were amorphous yet structurally homo- 
geneous; the thickness ranges from 15 to 40 nm. DC transport measurements were used to 
characterize T c . 

For the samples with a screening metal plane (Fig. H]) we used the following procedure: A 
30nm thick Al film was deposited on an insulating substrate. The Al was allowed to oxidize 
in ambient conditions for 24 hours. This process has been shown to produce ~ 5nm thick 



pinhole free high quality AlO insulators 21J. The T c of the Al layer was 1.1K, well below 
the measurement temperature. A disordered SC film was then deposited on the AlO film. 
As a cross-check, three additional pairs of "screened" and "bare" samples were produced 
on different substrates (SiO, sapphire or MgO) all showing similar differences in T c . In 
particular a pair of samples evaporated on the same substrate, Si + Ag + AlO, having 
different thickness of the AlO insulator exhibited a lower Tc for the sample having a thicker 
insulator. 

THz spectroscopy was used in the past to confirm the BCS theory since it probes the 



energy range of the superconducting gap [3514371] . Optical methods offer a number of crucial 



advantages over DC and tunneling measurements. In the first place a macroscopic area 
is irradiated homogeneously through the whole sample; thus it is a true volume averaging 
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measurement unlike DC methods that are more sensitive to actual current trajectories and 
film morphology. Secondly it relinquishes the need for metal contacts which may affect the 
electron system by screening. Finally, THz spectroscopy can be employed to study also 

:>hase of the SIT, which are not accessible by tunneling 



samples very deep in the insulating 
methods. The experimental setup J40| is based on backward wave oscillators as powerful 
radiation sources to emit continuous- wave, coherent radiation which can be tuned over 
the frequency range of 0.05 — 1.1 THz, corresponding a photon energy of 0.18 — 4.5 meV. 
The beam traversing through the sample is subject to multi-reflection on the two parallel 
boundaries of the substrate (i.e. Fabry-Perot interference) which strongly increases the 
interaction with the film. Subsequently it is detected by a Golay cell or He-cooled bolometer. 
A reference arm completes the arrangement for a Mach-Zehnder interferometer in order to 
measure the phase shift in addition to the sample transmission for each frequency. The 
specimens were mounted in an optical 4 He bath cryostat, allowing to cool down to T = 
1.85 K. 

For each temperature and applied frequency we measured simultaneously the transmis- 
sion, Tr, and phase shift, 0, of the radiation passing through the studied sample, i.e. the 
film on the substrate; the bare substrate was characterized beforehand. Obtaining these 
two independent parameters enabled us to directly extract the complex response function 
(e.g. the real and imaginary parts of the complex conductivity, o\ and o"2, respectively) as 
demonstrated in Fig. [5] [39( . 

There are two ways to determine the energy gap A from the raw data. In the first 
approach, Tr(f) and (j)(f) are fitted by Fresnel's equations for a thin film on a transparent 



substrate |3§J with frequency dependent optical parameters derived from microscopic models. 
In the normal state the maximum of the Fabry-Perot oscillations hardly vary with frequency; 
the conductivity is given by the Drude model with the DC conductivity and the scattering 
rate as fitting parameters. For T < T c both Tr and <fi change dramatically with /. In 
this case the THz data are fitted to the Mattis-Bardeen equations complemented by a finite 
scattering rate [32]; the gap 2 A and the plasma frequency uj p are the fitting parameters. 
This procedure is demonstrated in Fig. [5^ for a film of NbN. 

The second analysis procedure does not rely on any microscopic model. From the two 
independent quantities Tr and (f) the complex conductivity, o\ and 02, is calculated via the 



Fresnel equations (for multiple reflections) [38J . From the frequency dependent conductivity, 
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FIG. 5: Obtaining the energy gap value and the complex conductivity from the trans- 
mission and phase-shift data, a, <f>(f)c/f and Transmission versus energy (frequency) of a NbN 
film having a critical temperature of 15 K. The thick curves are the raw data and the thin curves 
are the fits to Mattis-Bardeen. The growing deviation of the phase shift over frequency from the 
normal state at T < T c with lowering T or / indicates the increasing superfluid stiffness. The peak 
in the transmission is a precursor for the value of 2 A. b, Real part of the complex conductivity 
{a{) as a function of temperature and energy, c, Imaginary part of the complex conductivity (02) 
as a function of temperature and energy. 

the value of 2 A can be easily extracted by the kink in &i(f). Photon energies hf > 2 A 
allow excitations across the superconducting gap, leading to a sharp increase in absorption. 
Radiation with smaller frequency (hf < 2A) can not break Cooper pairs, and only normal 
quasi-particle absorption is possible. With lowering the temperature, quasi-particle excita- 
tions become sparse and <J\(f) remains flat down to lowest frequencies or increase with a 
further decrease in frequency - depending on the value of the scattering rate. a 2 (which 
mainly depends on </>//) diverges as / — > 0; according to the Kramers-Kronig relation, this 
corresponds to a strong zero frequency peak, representing superconductivity. This analysis 
procedure was employed to produce the plots in Fig. [5] panels b and c. 
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